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ABSTRACT 

We present sensitive molecular line observations of the metal-poor blue compact dwarf I Zw 18 
obtained with the IRAM Plateau dc Burc interferometer. These data constrain the CO J = 1 — » 
luminosity within our 300 pc (FWHM) beam to be Leo < 1 x 10 5 K km s _1 pc 2 (Ico < 1 K km s -1 ), an 
order of magnitude lower than previous limits. Although I Zw 18 is starbursting, it has a CO luminosity 
similar to or less than nearby low-mass irregulars (e.g. NGC 1569, the SMC, and NGC 6822). There is 
less CO in I Zw 18 relative to its B-band luminosity, H I mass, or star formation rate than in spiral or 
dwarf starburst galaxies (including the nearby dwarf starburst IC 10). Comparing the star formation 
rate to our CO upper limit reveals that unless molecular gas forms stars much more efficiently in 
I Zw 18 than in our own galaxy, it must have a very low CO-to-H2 ratio, ~ 10 -2 times the Galactic 
value. We detect 3mm continuum emission, presumably due to thermal dust and free-free emission, 
towards the radio peak. 

Subject headings: galaxies: individual (I Zw 18); galaxies: ISM; galaxies: dwarf, radio lines: ISM 
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1. INTRODUCTION 

With the lowest nebular metallicity in the nearby un i- 
verse (12 + log O/H w 7.2. lSkillman k Kennicut^fl993l) . 
the blue compact dwarf I Zw 18 plays an important role 
in our understanding of galaxy evolution. Vigorous ongo- 
ing star formation implies the prese nce of molecular gas, 
but di rect evidence has been elusive. IVidal-Madiar et"aLl 
(2000) showed that t here is not significant diffuse H 2 , but 
ICannon et"aH (|2002f ) found ~ 10 3 M w of dust organized 
in clu mps with sizes 50 - 100 pc. IVidal-Madiar et all 
(|2000l) did not rule out co mpact, dense molecular clouds, 
and lCannon et all (|2002l ) argued that this dust may in- 
dicate the presence of m olecular gas. 

Observations by Arnault et al.1 (|1988f ) and 
iGondhalekar et all dl99J) failed to detect CO J = 1 -» 
emission, the most commonly used tracer of H2. This 
is not surprising. The low dust abundance and intense 
radiation fields found in I Zw 18 may have a dramatic 
impact on the formation of H2 and structure of molecular 
clouds. A large fraction of the H2 may exist in extended 
envelopes surrounding relatively compact cold cores. In 
these envelopes, H? self- shields while CO is dissociated 
(jMalonev k Black! Tl988l ). The result may be that in 
such galaxies [CII] or FIR emission trace H2 better tha n 
CO dMadden et a.lJfl997t ITsrael Il997at [Pak et al.l[l99cl . 
Further, H2 may simply be underabundant, as there is a 
lack of grains on which to form while photod issociation 
is enh anced by an intense UV field. Indeed, iBell et aTl 
(|2006[ ) found that at Z = Z®/W0, a molecular cloud 
may take as long as a Gyr to reach chemical equilibrium. 

A low CO content in I Zw 18 is then expected, and a 
stringent upper limit would lend observational support to 
predictions for molecular cloud structure at low metallic- 
ity. However, while the existing upper limits are sensitive 
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in an absolute sense, they do not even show I Zw 18 to 
have a lower normalized CO content than a spiral galaxy 
(e.g. less CO per £?-band luminosity). Th e low luminos- 
ity (M B w -14.7. iGil de Paz et alTl2003j) and large dis- 
tance (d=14 Mpc. llzotov k Thuanl [2004) of this system 
require very sensitive observations to set a meaningful 
upper limit. 

In this letter we present observations, obtained with 
the IRAM Plateau de Bure Interferometer (PdBI) 5 , that 
constrain the CO luminosity, LcO: to be equal to or less 
than that of nearby CO-poor (non-starbursting) dwarf 
irregulars. 

2. OBSERVATIONS 

I Zw 18 was observed with the IRAM Plateau de 
Bure Interferometer on 17, 21, and 27 April and 13 
May 2004 for a total of 11 hours. The phase calibrators 
were 0836+710 (F„(115GHz) « 1.1 Jy), and 0954+556 
(^(llSGHz) « 0.35 Jy). One or more calibrators with 
known fluxes were also observed during each track. The 
data were reduced at the IRAM facility in Grenoble us- 
ing the GILDAS software package; maps were prepared 
using AIPS. The final CO J = 1 — > data cube has beam 
size 5.59" x 3.42", and a velocity (frequency) resolution 
of 6.5 km s _1 (2.5 MHz). The velocity coverage stretches 
from vlsr ~ 50 to 1450 km s _1 . The data have an RMS 
noise of 3.77 mJy beam -1 (18 mK; 1 Jy beam -1 = 4.8 
K). The 44" (FWHM) primary beam completely covers 
the galaxy. Based on variation of the relative fluxes of 
the calibrators, we estimate the gain uncertainty to be 
< 15%. 

3. RESULTS 

3.1. Upper Limit on CO Emission 

To search for significant CO emission, we smooth the 
cube to 20 km s _1 velocity resolution, a typical line widt h 
for CO at our spatial resolution (e.g.. [Heifer et al1 [2003). 

5 Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain)." 
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Fig. 1. — CO 1 — > spectra of I Zw 18 towards the ra- 
dio continuum/Ha peak (left) and the highest significance spectra 
(right), which is still too faint to classify as more than marginal. 
The locations of both spectra are shown in Figure [2] Dashed hor- 
izontal lines show the magnitude of the RMS noise. 



The noise per channel map in this smoothed cube is 
(J20 ~ 0.25 K km s" 1 . Over the H I velocity range (710 
- 810 km s~ 1 . lvan Zee et aT1ll998f ). there are no regions 
with Ico,20 > 1 K km s _1 (4er) within the primary beam. 
We pick a slightly conservative upper limit for two rea- 
sons. First, if there were CO emission with this intensity 
we would be certain of detecting it. Second, the noise in 
the cube is slightly non-Gaussian, so that the false posi- 
tive rate for ico, 20 > 1 K km s _1 — estimated from the 
negatives and the channel maps outside the H I velocity 
range — is ~ 0.2%, very close to that of a 3a deviate. 

For d = 14 Mpc, the synthesized beam has a FWHM 
of 300 pc and an area of 1.0 x 10 5 pc 2 . Our intensity 
limit, Ico < 1 K km s , therefore translates to a CO 
luminosity limit of Leo < 1 x 10 5 K km s _1 pc 2 . 

There is a marginal signal toward the southern knot 
of Ha emission (9 h 34 m 02 s .4, 55°14'23".0). This emission 
has the largest | ico. 20 1 found over the H I velocity range, 
corresponding to Leo ~ 8x 10 4 K km s" 1 pc 2 , just below 
our limit. This same line of sight also shows \Ico\ > 2cr 
over three consecutive channels, a feature seen along only 
one other line of sight (in negative) over the H I velocity 
range. The marginal signal is sugges tively l ocate d in the 
southeast of I Zw 18, where [Cannon et al.l (|2002f ) identi- 
fied several potential sites of molecular gas from regions 
of relatively high extinction. While tantalizing, the sig- 
nal is not strong enough to be categorized as a detection. 
Figure [Q sh ows CO spectra towards the Ha/radio contin- 
uum peak (jCannon et al.l l2002. 2005. lHunt et al.ll2005aL 
see Figure [5]) and this marginal signal. 

3.2. Continuum Emission 

We average the data over all channels and produce a 
continuum map with noise ctiisghz = 0.35 mJy beam -1 . 
The highest value in the map is Llisghz = 1.06±0.35 mJy 



beam" 1 at a 200 o = 9 h 34 m 02 s .l, <5 2 ooo = +55° 14' 27".0. 
This is withi n a fraction of a beam of the 1.4 GHz peak 
identified bv lCannon et~afl (120051 a.nnn = 9 h 3 4 m 02M, 
£2000 = +55° 14' 28".06) and lHunt et all (|2005al a 20 oo = 
9 h 34 m 02 s , <5 2 ooo = +55° 14' 29".06). Figure [2] shows the 
radio continuum peak and 115 GHz conti nuum contours 
plott ed over Ha emission from I Zw 18 (| Cannon et al.l 
2002). There is only one other region with |/n5GHz| > 
3(Tii5GHz within the primary beam and the star-forming 
extent of I Zw 18 occupies « 10 % of the primary beam. 
Therefore, we estimate the chance of a false positive co- 
incident with the galaxy to be only ~ 10%. 



4. DISCUSSION 

Here we discuss the implications of our CO upper limit 
and continuum detection. We adopt the following prop- 
erties f or I Zw 18, all scaled to d = 14 Mpc: Mb = 
-14.7 dGil de Paz et al.l l20(jl. M HI = 1.4 x 10 8 M 
(Ivan Zee et al. 19981) . Ha luminosity logqn H a = 39.9 erg 
s" 1 ([Cannon et al.l2002tlGil de Paz e t al. 2003h . 1.4 GHz 
flux F 1A = 1.79 mJy (jCannon et al.ll2005h . 



4.1. Point Source Luminosity 

Our upper limit along each line of sight, Leo < 
1 x 10 5 K km s _1 pc 2 , matches the luminosi ty of 
a fai rly massive Galactic giant molecular cloud (jBlitzl 
1993). For a Galactic CO-to-H 2 conversion factor, 2 x 
10 20 cm~ 2 (K km s -1 ) -1 , the corresponding molecular 
gas mass is MmoI ~ 4.4 x 10 5 M©, s imilar to the mass of 
the Orion-Monoceros complex (e.g. IWilson et aLH "2005). 

4.2. Comparison With More Luminous Galaxies 

In galaxies detected by CO surveys, the CO content 
per unit B-band luminosity is fairly constant. Figure [3] 
shows the CO luminosity normalized by B-band lumi- 
nosity, Loo I Lb, as a function of absolute B-band mag- 
nitude {Lb is extinction corrected). Leo /Lb is nearly 
constant over two orders of magnitude in Lb, though 
with substantial scatter (much of it due to the extrapo- 
lation from a single pointing to Leo)- 

Based on these data and assuming that Leo is not 
a function of the metallicity of the galaxy, we may ex- 
trapolate to an expected CO luminosity for I Zw 18. 
For Afs,iZwi8 ~ —14.7 the CO luminosity correspond- 
ing to the median value of Leo/ Lb (dashed line) in 
Figure [3] is L C o,iZwl8 ~ 1-7 x 10 6 K km s" 1 pc 2 . 
The Ha, 1.4 GHz , and H I lumin osities lead to simi- 
lar predictions. lYoung et al.1 (|1996| ) found M# 2 /L# Q rs 
10L©/M© for Sd-Irr gala xies, wh ich impli es L cojzwis ~ 
4 x 10 6 K km s _1 pc 2 . iMurgia et~all (|2005h measured 
Fgo/Fia ~ 10 Jykms -1 (mJy) -1 for spirals, that would 
imply Lco,iZwi8 ~ 10 7 K km s" 1 . For Sd/S m galax- 
ies, M H2 /M HI « 0.2 (|Young fc Scovillelll991h . leading 
to L co ,iZwi8 ~ 5 x 10 6 K km s _1 pc 2 . Both M H2 /L Ha 
and Mh2/Mhi tend to be even higher in earlier-type 
spirals. 

Therefore, surveys would predict Lcojzwis ^ 2 x 
10 6 K km s _1 pc 2 , very close to the previou sly established 
upper limits of 2 — 3 x 10 6 K km s _1 pc 2 ([Arnault et al.1 
ll988tlGondhalekar et aDll998f ). With the present obser- 
vations, we constrain Leo < 1 X 10 5 K km s _1 pc 2 and 
thus clearly rule out Leo ~ 10 6 K km s _1 pc 2 . This 
may be seen in Figure 02 even if I Zw 18 has the highest 
possible CO content, it will still have a lower Leo/ Lb 
than 97% of the survey galaxies. 

4.3. Comparison With Nearby Metal-Poor Dwarfs 

The subset of irregular galaxies detected by CO surveys 
tend to be CO-rich and actively star-formin g, resembling 
scaled-down versions of spir al galaxies ([Young et al.1 
119951 119961 : iLerov et all 120051 ). Such galaxies may not 
be representative of all dwarfs. Because they are nearby, 
several of the closest dwarf irregulars have been detected 
despite very small Leo- With their low masses and 
metallicities, they may represent good points of compar- 
ison for I Zw 18. Table [T] and Figure [3] show CO lumi- 
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Fig. 2.— V-band and Ha (right, Cannon ct al. 2002) images of I Zw 18. Overlays on the left image show the size of the synthesized 
beam and the locations of the spectra shown in Figure [T] Contours on the right image show continuum emission in increments of 0.5<r 
significance and the location of the radio continuum peak. The primary beam is larger than the area shown. Both optical maps are on 
linear stretches. V-band data obtained from the MAST Archive, originally observed for GO program 9400, PI: T. Thuan). 



nosities and Leo/ Lb for four nearby dwarfs: NGC 1569, 
the Small Magellanic Cloud (SMC), NGC 6822, and 
IC 10. The SMC, NGC 1569, and NGC 6822 have 
Loo ~ 10 5 K km s _1 pc 2 , close to our upper limit, and 
occupy a region of Leo / Lb-Lb parameter space similar 
to I Zw 18. All four of these galaxies have active star for- 
mation but very low CO content relative to their other 
properties. 

We test whether our observations would have detected 
CO in NGC 1569, the SMC, and IC 10 at the plausible 
lower limit of 10 Mpc (from Hq = 72 km s _1 ) or our 
adopted distance of 14 Mpc. We convolve the integrated 
intensity maps to resolutions of 210 and 300 pc and mea- 
sure the peak integrated intensity. The results appear in 
columns 4 and 5 of Table [T] The PdBI observations of 
NGC 1569 resolve out most of the flux, so we also apply 
this test t o a distribution wit h the size and luminosity 
derived by iGreve et al.l (|1996h from single dish observa- 
tions. Our observations would detect an analog to IC 10 
but not the SMC, with NGC 1569 an intermediate case. 
With a factor of ~ 3 better sensitivity (requiring ~ 10 
times more observing time) we would expect to detect 
all three nearby galaxies. However, achieving such sen- 
sitivity with present instrumentation will be quite chal- 
lenging. ALMA will likely be necessary to place stronger 
constraints on CO in galaxies like I Zw 18. 

IC 10 may be the nearest blue compact dwarf 
(|Richer et al.ll200"l[ ). so it may be telling that we would 
detect it at the distance of I Zw 18. The blue compact 
galaxies that ha ve been detected in CO have Leo /Lb 
similar to IC 10 (|Gondhalekar et aL| [l998. the diamonds 
in Figure [3]). Most searches for CO towards BCDs have 
yielded nondetections, so those detected may not be rep- 
resentative, but I Zw 18 is clearly not among the "CO- 
rich" portion of the BCD population. 



4.4. Interpretation of the Continuum 



We measure continuum intensity of i r ii5GHz = 1-06 ± 
0.35 mJy towards the radio continuum peak. The 
continuum is detected along only one line of sight, 
so we refer to it here as a point source and com- 
pare it to integrated values for I Zw 18. -FiisgHz is 
expected to be the product of mainly two types of 
emission: thermal free-free emission and thermal dust 
emission. At long wavelengths, the integrated ther- 
mal free- f ree emission is fi 4GHz (free — free) » 0.52 - 
0.75 mJy (|Cannon et alj|2005t [Hunt et al.ll2005a|) . imply- 
ing i 7 'ii5GHz(free — free) = 0.36 - 0.51 mJy at 115 GHz 
[F v ex v~ a l ). The Ha flux pred icts a similar value , 
^115GHz (free -free) = . 34 mJ y (jCannon et all 120051 . 
Equation 1). iHunt et all ((2005b) placed an upper limit 
of F„(850) < 2.5 mJy on dust continuum emission at 
850/im; t his is consistent with the ~ 5 x 10 3 M Q esti- 
mated by I Cannon et alj (|2002f ) given almost any reason- 
able dust properties. Extrapolating this to 2.6 mm as- 
suming a pure blackbody spectrum, the shallowest plau- 
sible SED, constrains thermal emission from dust to be 
< 0.25 mJy at 115 GHz. Based on these data, we would 
predict FhsgHz ^5 0-75 mJy. Thus our measured Fh 5 ghz 
is consistent with, but somewhat higher than, the ther- 
mal free-free plus dust emission expected based on opti- 
cal, centimeter, and submillimeter data. 

4.5. Relation to Star Formation 

I Zw 18 has a star formation rate ~ 0.06 - 0.1 M Q 
yr _1 , based on Ha and cm radio continuum measure- 
ments ([Cannon et al.|[2002l : lKennicutt 1998a; Hun t et aLl 
l2005al) . Our continuum flux suggest s a slightly higher - 
value k, 0.15 - 0.2 M yr" 1 (following IHunt et alJl2Q05a[ 
lCondon|[l99 2), with the exact value depending on the 
contribution from thermal dust emission. For any value 
in this range, the star formation rate per CO luminosity, 
SFR/Lqo is much higher in I Zw 18 than in spirals. For 
comparison, our upper limit and the molecular "Schmidt 
Law" derived bv lMurgia et a l. (2002) predicts a star for- 



4 



10° 



J 10 " 



10" 



10" 



10"" 



10" 



UM46: 
U*45l 



• • !t ••• . 



NGC5251* 



Market; o». Hr°5> 

UM48& • * 



IZwlB 



X Nearby low mass dwarfs 
O Blue compact galaxies 
(Gondhalekar et al. 1998) 
I . . . I . . . I . . . I 



-14 



-16 



-18 -20 



-22 



-24 



Fig. 3. — CO luminosity normalized by absolute blue mag- 
nitud e for galaxies with Hubble Type Sb o r later (black cir- 
cles. lYoung et all IT995I ; lElfhag et aLl Il99l IBoker et all [20031 ; 
ILerov et al.l 120051) . We also plot nearby dwarfs from Ta- 
ble [T] (crosses) and blue compact galaxies compiled by 
I Gondha lekar ct al. (1998, , diamonds). The shaded regions shows 
our upper limit for I Zw 18, with the range in Mb for distances from 
10 to 20 Mpc. The dashed line and light shaded region show the 
median value and la scatter in Leo/ Lb f° r spirals and dwarf star- 
bursts. Methodology: We extrapolate from Iqq in central pointings 
to Lqq assumi ng the CO to have an exponential profile with scale 
length 0.1 d,2s l|Young ct al. 1995), including only galaxies where 
the central pointing measures > 20% of Leo- We adopt B mag- 
nitudes (corrected for internal and Galactic extinction), distances 
(Tully-Fisher when available, otherwi se Virgocentric-flow corrected 
Hubble flow), and radii from LEDA UPaturel et al.ll2003T) . 



matio n rat e < 2 x 10 ~ 4 Mm yr" 1 . Fits by lYoung et al.l 
(1996) and iKennicuttj (|1998bl applied to just the molec- 
ular limit) yield similar values. Again, I Zw 18 is similar 
to the SMC and NGC 6822, which have st ar formation 
rates of 0.05 M yr" 1 and 0.04 M yr" 1 (|Wilke et all 
[200l Iisraellfl997bh and L co ~ 10 5 K km s -1 pc 2 . 



4.6. Variations in Xco 

Several calibrations of the CO-to-H2 conversion factor, 
Xqo as a function of metallicity exist in the literature. 
The topic has been controversial and thes e calibrations 
range from little or no dependence (e.g. iWalterl [20031 : 
iRosolowskv et al J 12003ft to v ery steep dependence (e.g., 
Xco oc Z~ 2 7 llsraell Il997aj) . Comparing the star for- 
mation rate to our CO upper limit, we may rule out 



that I Zw 18 has a Galactic Xco unless molecular gas 
in I Zw 18 forms stars much more efficiently than in the 
Galaxy. Either the ratio of CO-to-H2 is low in I Zw 18 
or molecular gas in this galaxy forms stars with an effi- 
ciency two orders of magnitude higher than that in spiral 
galaxies. 

5. CONCLUSIONS 

We present new, sensitive observations of the metal- 
poor dwarf galaxy I Zw 18 at 3 mm using the Plateau de 
Bure Interferometer. These data constrain the integrated 
CO J = 1 — > intensity to be lco < IK km s _1 over our 
300 pc (FWHM) beam and the luminosity to be Lco < 
1 x 10 5 K km s" 1 pc 2 . 

I Zw 18 has less CO relative to its B-band luminosity, 
H I mass, or SFR than spiral galaxies or dwarf starbursts, 
including more metal-ri ch blue compact galaxies such as 
IC 10 (Z IC io ~ Zq/A, iLeTeTall [2003ft . Because of its 
small size and large distance, these are the first observa- 
tions to impose this constraint. 

We show that I Zw 18 should be grouped with several 
local analogs — NGC 1569, the SMC, NGC 6822 — as 
a galaxy with active star formation but a very low CO 
content relative to its other properties. In these galax- 
ies, observations suggest that the environment affects the 
molecular gas and these data suggest that the same is 
true in I Zw 18. A simple comparison of star formation 
rate to CO content shows that this must be true at a 
basic level: either the ratio of CO to H2 is dramatically 
low in I Zw 18 or molecular gas in this galaxy forms stars 
with an efficiency two orders of magnitude higher than 
that in spiral galaxies. 

We detect 3mm continuum with Fn 5 GHz = 1-06 ± 
0.35 mjy coincident wit h the radio peak id entified by 
ICannon et all (|2005ft and lHunt et~aT1 (|2005aft . This flux 
is consistent with but somewhat higher than the thermal 
free-free plus dust emission one would predict based on 
centimeter, submillimeter, and optical measurements. 

Finally, we note that improving on this limit with cur- 
rent instrumentation will be quite challenging. The order 
of magnitude increase in sensitivity from ALMA will be 
needed to place stronger constraints on CO in galaxies 
like I Zw 18. 
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TABLE 1 
CO in Nearby Low Mass Galaxies 



Galaxy 


M B 


^CO 


ICO, 210 a 


ICO, 300 a 


Reference 




(mag) 


(K km s -1 pc 2 ) 


(K km s" 1 ) 


(K km s- 1 ) 




NGC 1569 


-16.5 


1.2 X 10 5 


1.1 


0.8 


Greve et al. (1996) 




-16.5 


0.2 X 10 s 


0.8 


0.5 


Tavlor et al. ( 1999) 


SMC 


-16 


1.5 x 10 5 


0.5 


0.1 


Mizuno et al. (2001. 2006") 


NGC 6822 


-16 


1.2 x 10 5 






Israel (1997b) 


IC 10 


-16.5 


2.2 x 10 6 


3.8 


2.2 


Lcroy et al. (2006) 


I Zw 18 


-14.7 


< 2 x 10" 






Arnault et al. (1988); Gondhalckar et al. (1998) 


I Zw 18 


-14.7 


< 1 x 10 5 


< 1 


< 1 


this paper 



a Peak integrated intensity at 210 and 300 pc, corresponding to our beam 
size at 10 and 14 Mpc, respectively. 
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